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A NUMERICAL METHOD FOR THE STRESS ANALYSIS OF STIFFENED-SHELL STRUCTURES 

UNDER NONUNIFORM TEMPERATURE DISTRIBUTIONS 1 

By Rickard R. Heldexfbds 


SUMMARY 

A numerical method is presented for the stress analysis 
of stiff ened-shell structures of arbitrary cross section under 
nonuniform temperature distributions. The method is based 
on a previously published procedure that is extended to include 
temperature effects and multicell construction. The applica- 
tion of the method to practical problems is discussed, and an 
illustrative analysis is presented of a two-cell box beam under 
the combined action of vertical loads and a nonuniform tem- 
perature distribution. 

INTRODUCTION 

The effects of nonuniform temperature distributions, 
such as those produced by aerodynamic heating, are becoming 
of greater concern in the design of modem high-speed aircraft. 
The structural effects of temperature changes and the 
results of some analyses of a simplified, structure under 
nonuniform distributions of temperature have been dis- 
cussed in reference 1. The analytical methods considered 
in reference 1 were found, however, to yield inaccurate values 
for the secondary stresses in complicated structures, and 
in such cases some type of numerical approach is desirable. 
Numerical methods, however, usually require extensive and 
tedious calculations and they should be used only when 
satisfactory results cannot be obtained from a simplified 
analysis. 

Several numerical methods of stress analysis have been 
presented in the literature, but none contains provisions for 
temperature changes. In the present report, one. such 
method, the numerical procedure of reference 2, has been 
extended to include the effects of a honuniform distribution 
of temperature. In addition, the equations developed 
permit the analysis of a stiffened-shell structure of arbitrary 
cross section with any number of internal cells. The appli- 
cation of the method is discussed and illustrated by analysis 
of a two-cell box beam under the combined action of vertical 
loads and a nonuniform temperature distribution. 

DESCRIPTION OF THE NUMERICAL METHOD 
BASIC THEOBT 

The structure analyzed is an idealized representation of a 
multicell stiffened-shell structure (see fig. 1) and has the 


following characteristics: 

(1) The basic unit is a rectangular panel bounded on two 
parallel sides by extensionally flexible stringers and on the 
other two sides by rigid bulkheads. 

(2) The panels consist of sheet material and are assumed 
to carry shear stress only. The shear stress is constant with- 
in a given panel. 

(3) The stringers run parallel to the direction of the pri- 
mary stresses and carry axial load only. 

(4) The bulkheads lie perpendicular to the stringers and 
are rigid in their own plane but offer no resistance to warping 
out of their plane. 

(5) The structure i? loaded only at the bulkheads. 

(6) Material properties, cross-sectional dimensions, and 
temperature distribution do not vary along the length of a 
given bay. • 

With these assumptions about the basic elements of the 
structure, any type of stiffened shell can be analyzed, pro- 
vided taper is excluded. The state of stress in such a struc- 
ture can then be defined by suitable stress-strain relations 
and two types of displacements: 

(1} Stringer displacements, which are displacements, at 
the end of a bay, of each flexible stringer in a direction parallel 
to the stringer 

(2) Bay displacements, which are translations and rota- 
tions of the plane of each cross section defined by the rigid 
bulkheads 



1 Supersedes NACA TN 2241, “A Numerical Method for the Stress Analysis of Stfflened-SheU Structures under Nonunlform Temperature Efitrtbutfcms” hy Efchard B. Heldenfels, 1640. 
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Once the stress-strain relations are established for the. 
components of the idealized structure, equations of equilib- 
rium can be used to obtain relationships between the dis- 
placements. The equilibrium equation for the forces on an 
individual stringer yields an expression for the stringer dis- 
placement of any panel point in terms of the surrounding 
stringer displacements and the displacements of the two 
adjacent bays. From this general expression, equations 
equal in number to the unknown stringer displacements are 
obtained. The additional equations required for the deter- 
mination of the bay displacements are obtained from the 
equations of equilibrium of the shear forces on the cross sec- 
tions. These equations then completely define the dis- 
placements of the structure. In most cases the number of 
equations is so large that a direct solution would be imprac- 
tical and it has been found expedient to solve them by the 
recommended iteration procedure described in the next sec- 
tion. The required equations are derived in detail in the 
appendix. 

SOLUTION OF EQUATIONS BY ITERATION 

Matrix iteration often provides the easiest and quickest 
solution to the equations, and the procedure recommended 
is as follows: . . 

The equations to be solved can be written in matrix 
notation as 

[E\{d} = {c} (1) 

For purposes of iteration, these equations are rearranged to 
give 

{d}=[C]{d} + {6} (2) 

where 

[C] = [?7]-[5j . . ._ 

[. B ] square matrix of coefficients of general equations with 

diagonal terms reduced to unity 
\U] unit matrix 

{ d } column matrix of stringer and bay displacements 
{c} column matrix of constant terms in general equation; 
these terms arise from applied load and thermal 
expansion 

Initial approximate values of stringer and bay displace- 
ments {d 0 } are then selected. These. values may be deter- 
mined in any convenient manner; however, subsequent 
operations can be simplified, as explained in the appendix, if 
these values are chosen to correspond to elementary theory. 
Next, the initial displacement values are substituted into the 
righthand side of equation (2) to obtain a second approxi- 
mation { di } to the displacements 

{4M<71{d,} + {c} (3) 

and the differences between the second and initial approxi- 
mate displacement values are computed from the equation 

{ Adi } = { di } — {do} (4) 


The iteration process is then begun by using those dis- 
placement differences. The nth difference is defined as 

{Ad„} = {d n } — {d 0 } (5) 

and it can be easily verified that the use of those differences 
leads to the following matrix equation: 

{Ad« } = [(7]{Ad n _i} + {Adi} (6) 

The iteration process consists of a series of solutions of equa- 
tion (6), each successive solution yielding a better approxi- 
mation to the displacement differences than the previous 
one. The process is continued until successive solutions of 
equation (6) yield the same result, that is, until 

{Ad,} = {Ad,_ 1 } (7) 

The final displacements are then determined from the final 
differences by using equation (5) and the initial values. 

When equation (6) is being iterated, improved values 
should be used as soon as they are obtained; that is, each 
individual difference Ad, should be substituted into the 
{ Ad„_i } matrix immediately aft or calculation rather than 
at the end of the cycle. In this manner, each new value 
determined receives the benefit of all previous work and 
convergence is speeded. 

The iteration of differences reduces the work required to 
obtain a solution because smaller numbers are involved. 
However, it is essential that no errors be made in the. deter- 
mination of the first differences { Adi } since a singlo significant 
error will render the whole solution useless. 

CONVERGENCE OF THE ITERATION PROCESS 

In order to obtain more rapid convergence of the iteration 
process, - bay displacements and loads arc referred to the 
principal shear axes of each bay. The use of these axes 
greatly simplifies the equations for bay displacements by 
making each bay displacement independent of all other bay 
displacements and thus a function of the stringer displace- 
ments only. In addition, a special correction cycle is 
periodically introduced to bring the stringer forces on each 
cross section into equilibrium with the applied loads. Mathe- 
matically, the correction cycle is a special cycle that usos a 
certain combination of the basic equations. Its success in 
the particular case of the numerical method of stress analysis 
is a result of its physical significance, and in that respect it is 
similar to Southwell’s “group relaxations” (reference 3). 

The optimum frequency of application of the correction 
cycle depends largely on the characteristics of each individual 
problem and must be determined on a .basis of experience 
with the method. If this frequency cannot be determined 
from previous experience, it can bo approximated satisfac- 
torily by one that permits the disturbances to spread their 
significant effect over the structure between correction cycles. 



STRESS ANALYSIS OF STIFFEN ED -SHELL STRUCTURES 


1021 


The application of the correction cycle begins at a station 
where the displacements are known and then proceeds out- 
board. The corrections required to bring the first bay into 
equifibrium are determined, and the stringer displacements 
at its outboard end are changed accordingly before the cor- 
rections required by the second bay are calculated. 

EFFECT OF INTRODUCING NONUNTFOBM TEMPERATURE DISTRIBUTIONS 

The preceding method is applicable to any type of stress 
problem. iSTon uniform temperature distributions do not 
affect the general procedure but merely change the details. 
These effects are of two types: A change in the effective 
structure due to changes in elastic properties of the material 
with temperature and thermal stresses resulting from re- 
strained thermal expansion. The changes in elastic prop- 
erties are easily handled if the moduli are treated as variables 
during the derivation of the equations. Their effect is 
analogous to that of variations in stringer area and panel 
thickness. The presence of thermal expansion requires 
modification of the stress-strain relationships for the stringers 
but does not affect those for the panels. The equations for 
stringer displacements contain thermal-expansion terms that 
are analogous to the applied-load terms. Bay-displacement 
equations are unaffected by thermal expansion, but thermal- 
expansion tpxms appear in the equations used for the 
correction cycle. If a difference solution is iterated, the 
elementary solution should include the distributions of 
thermal strain associated with the primary thermal stresses, 
which may be obtained from the equations derived in the 
appendix: 

DISCUSSION OF THE NUMERICAL METHOD 

The application of a method, such as that just described, 
always poses a number of questions; for example, what are 
some of the limitations of the method, would it be advan- 
tageous to use some other method of analysis, and how should 
the structure be idealized? Some of the factors requiring 
consideration, other than those mentioned in the previous 
section, are therefore now discussed. 

VALIDITY OF BASIC ASSUMPTIONS 

The assumptions upon which the method is based are 
commonly accepted in the analysis of stiffened shells. Com- 
parison of theoretical and experimental results has estab- 
lished the fact that these assumptions will yield good results 
in most cases. Two important assumptions — that the 
bulkheads are rigid in their own plane and that the shear 
stress is constant in a given panel — may, however, introduce 
significant errors into the analysis in some cases. These 
assumptions are therefore examined in detail. 

The assumption of rigid bulkheads is satisfactory as long 
as the primary stresses run perpendicular to the bulkheads, 
hut, as demonstrated in reference 1, this assumption may 
not he good when dealing with problems involving thermal 


stress. Large temperature gradients along the length of the 
structure or across the depth of a bulkhead distort the real 
bulkhead and make the assumption of rigidity inapplicable. 
In many cases, however, these effects are small and the as- 
sumption yields satisfactory results. 

The numerical method could be extended to include the 
effects of bulkhead flexibility. Such an analysis, however, is 
very cumbersome and tedious and, if the equations are solved 
by iteration, the process is often very slowly convergent. 
Therefore, these extensions are not considered herein. 

The assumption of constant shear stress in a given panel 
simplifies the development of the equations, and it yields 
satisfactory results if the bulkheads are reasonably close 
together. Cases arise, however, in which the assumption 
will lead to unreasonable results because the assumed con- 
stant shear stress is a poor approximation to a shear stress 
which should he changing rapidly in the spanwise direction. 
This situation is usually accompanied by slow convergence 
of the iteration process. This difficulty, however, can be 
minimized by reducing the bulkhead spacing of the idealized 
structure since it occurs only when the total shear stiffness 
of the panels joined to a stringer exceeds the extensional 
stiffness of that stringer. 

IDEALIZATION' OF AN ACTUAL STRUCTURE 

The idealization process described in reference 2 is straight- 
forward. However, it provides an opportunity for the stress 
analyst to exercise his engineering judgment and thus to 
simplify the analysis. By restricting the analysis to only 
a part of the "structure or by using a comparatively simple 
idealized structure, tide time required for the analysis can he 
substantially reduced. Such simplifications, however, can 
reduce the value of the results, and a compromise between 
speed and exactness is required. 

The number and location of the idealized stringers com- 
pletely define the stress-distribution shapes obtainable from 
the analysis. (For example, in an idealized shell of n stringers, 
there are n possible types of independent normal-stress 
distributions, three of which can be determined from ele- 
mentary theory, the remaining n- 3 being statically indeter- 
minate.) Stringer location is thus an important part of the 
idealization process and in conventional problems the loca- 
tions should be selected after consideration of the charac- 
teristics of the actual structure, the nature of the expected 
results, and the time available for the analysis. When non- 
uniform temperature distributions are involved, the shape 
of the temperature distribution should also he considered 
because the thermal-stress and temperature distributions 
will have similar shapes and the analysis will yield good 
results only if the idealized structure permits a stress distri- 
bution of that shape. 

The bulkhead spacing usually is the same in the idealized 
and actual structures, but the idealized spacing should never 
be so .large that trouble is caused by the assumption of 
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constant shear stress. A proper bulkhead spacing is one for 
which the extensional stiffness of each stringer element is 
greater than the sum of the shear stiffnesses of the adjacent 
panels so that no negative terms appear on the right-hand 
side of the equations for the stringer displacements. 

CALCULATING PROCEDURES 

All the calculations required by the numerical procedure 
(determining the coefficients of the equations, solving by 
iteration, and computing the stresses) are routine and 
involve only simple arithmetic. The calculations can be 
easily arranged in tabular form so that the bulk of the work 
can be done by modem automatic computing machinery or 
by a computer who does not need to have a knowledge of the 
structural theories involved. 

In any problem that involves extensive numerical work, 
errors are very apt to occur. One of the advantages of the 
numerical method described herein is that a number of 
checking procedures can be devised to check the various 
steps in the calculations. No attempt is made to describe 
the many possible checks; a few, however, have been indi- 
cated in the illustrative example. 

Solution of the equations by simple iteration also possesses 
another advantage with regard to errors. Values obtained 
from successive cycles of iteration show trends that can 
be observed by an experienced computer, and errors can 
be detected by their effects on these trends. Errors that 
do appear during the iteration process eventually work 
themselves out but may adversely affect the rate of con- 
vergence. 

APPLICATION OF THE NUMERICAL METHOD 
DESCRIPTION OP THE PROBLEM 

The application of the numerical method is illustrated by 
an analysis of the idealized two-cell box beam shown in 
figure 2. The cross section is symmetrical about the hori- 
zontal center line and the beam is untapered; however, the 
stringer areas and sheet thicknesses vary from bay to bay. 

The box beam is loaded by four concentrated vertical loads 
applied at the bulkhead stations along the inner web; in addi- 
tion, it is subjected to the arbitrarily selected distribution 



Figure 2.— Idealised structure used In Illustrative example* (Stringer areas and skin thick- 
nesses are listed in tables T and II.) 


of temperature increase shown in figure 3. The temper- 
ature is highest at the tip and along the front web and 
decreases in both spanwise and chordwise directions, but it 
is constant across the depth of the beam. The beam is 
assumed to be constructed of 75S-T6 aluminum alloy which 
has the variation of elastic properties with temperature 
increase shown in figure 4. These data are the same as 
those used in reference 1 . 

It is assumed that no thermal stresses were present at 60° F. 
Since the method of analysis involves the assumption that 
no changes in temperature distribution occur over each 
element, the temperature used in the calculations was the 
temperature at the center of the element concerned. 




Fiqube 4. — Variation of elastic properties of 76S-T6 aluminum altor with temperature 
Increase (reference 1). 
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DETAILS OF THE ANALYSIS 

Since the structure and the temperature distribution are 
symmetrical about the horizontal center line, the analysis 
can be restricted to one cover. Two solutions are required 
however, to determine the total stress since the thermal- 
stress system is symmetrical about the horizontal center 
line, but the load-stress system is antisymmetries!. In'this 
way, the analysis requires the solution of two sets of equa- 
tions (one of 20, the other of 24) which can be solved more 
easily than the set of 44 needed for a single analysis of the 
complete box. 

The computations required are given in tabular form with 
most tables containing two parts — one related to the load 
stresses and the other related to the thermal stresses. The 
final solution is obtained by the superposition of these two 
solutions. The rectangular cross section and its symmetry 
permit several simplifications of the general equations. In 
each case the equations used are listed. The notation is 
described in the appendix. Methods used to check the 
calculations are also given in the tables. The checking 
methods used were determined from mathematical relation- 
ships existing between the coefficients of the equations and 
from equilibrium of forces. 

Tables I and IE present the physical characteristics and 
stiffness parameters of the individual stringers and panels. 
Table III gives the location of the principal shear axes of 
each bay and the coefficients of the bay-displacement equa- 
tions. The location of the principal inertia axes of each 
bay, the coefficients used in the correction cycle, andthe 
initial stringer displacements are given in table IV. The 
coefficients of the stringer-displacement equations are tabu- 
lated in table V. Table VT contains the [0] matrices used 
for the iteration. The rows and columns have been inter- 
changed in order that the matrix multiplications required 
will consist of the cumulative multiplication of the adjacent 
numbers in two columns. Table VIL is a similar arrange- 
ment of data required for the correction cycle and also lists 
each correction determined. The displacements obtained 
from each cycle of iteration are given in table VlLL and the 
correction cycles are indicated. Table IX contains the cal- 
culation of each type of stress and the superposition required 
to obtain the total stresses. 

The numerical calculations in this example were done by a 
computer who had previous experience with the method. 
The following times were required: 


Setting up the equations (table I to table VII) 3 days 

Solving the equations (table VIII) 4 days 

Computing stresses (table IX) 1 day 


In this example, the displacements were computed to six 
decimal places (five or six significant figures) in order that 
the stress would be accurate to 1 psi and thus would provide 
accurate equilibrium checks. Most practical problems will 
not require such numerical accuracy and a smaller number 
of decimal places should be used in order to speed the solu- 
tion. It is estimated that the time required to solve this 
example could have been reduced by one-half if the number 
of decimal, places had been reduced from six to four. This 
reduction would have given stresses accurate to 100 psi or 
about 1 percent of the maximum stress. 


RESULTS OF THE CALCULATION 

The results of the calculation are shown graphically in 
figure 5 by spanwise and chordwise plots of the stringer 
stresses in the top and bottom covers and a spanwise plot 
of shear stresses in the webs. The spanwise plots have a 
jagged appearance because stringer areSs and sheet thick- 
nesses •are assumed constant in each bay. with an abrupt 
change at the bulkheads. The dashed lines in the plots of 
stringer stresses are the values obtained from an elementary 
analysis. 

CONCLUDING REMARKS 

A numerical method for the stress analysis of stiff ened-shell 
structures under nonuniform, temperature distributions has 
been presented. The method is not applicable to the solu- 
tion of all structural problems involving temperature effects 
because it requires extensive and tedious calculations and 
because the basic assumptions of bulkheads rigid in their own 
plane and constant shear stress in a given panel occasionally 
Lead to unsatisfactory results. It is, however, a powerful 
tool for the solution of many structural problems because: 

(1) It is a means for accurately determining all types of 
secondary stresses in complicated structures that cannot be 
satisfactorily analyzed by simplified methods. 

(2) It is sufficiently flexible to cope with a wide variety of 
structural problems involving nonunifoim temperature dis- 
tributions. 

(3) It involves only simple arithmetic that can be handled 
by automatic computing machinery. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 12, 1960. 



(a) Spanwise distribution or upper-surface stringer stress. 
Figube 5. — Calculated stress distribution in the Idealized structure. 
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(b) Span wire distribution of lower-surface stringer stress. 
Ftoure 5.— Continued. 



(c) Cbordwlse distribution of upper-surface normal stress. 
Figure S.— Continued. 



(d) Cbordw lse distribution of lowebsorface normal stress. 
Figure S. — Continued. 



(e) Span wise distribution of web shear stress. 
Figure 6.— Concluded. 


APPENDIX 

DERIVATION OF GENERAL EQUATIONS 

The general equations required for the numerical analysis 
of a stiffened shell of arbitrary cross section with any number 
of internal cells and under a nonuniform temperaturcjdis- 
tribution are developed. The basic assumptions and a gen- 
eral description of the method have been given previously 
and are not repeated. 

SYMBOLS 

A cross-sectional area of stringer, square inches 

b width of panel on k grid line, inches 

E r ~- modulus of elasticity, psi 
F applied force, pounds 

G modulus of rigidity, psi 

h width of panel on grid line, inches 

I moment of inertia, inches* 

J shear stiffness parameters 

l, rrt coordinates of a special set of axes 
L ~ length of bay, inches 
M applied moment, inch-pounds 

P axial load in stringer, positive for tensile load, 

pounds 

Q area moment, inches® 

r normal distance to panel on k grid line, positive ill 

positive 2-direction, inches 

T temperature increment, measured from tempera- 

ture of zero thermal stress which is 60° F in the 
example presented, degrees Fahrenheit 
t panel thickness, inches 

u , v, w displacements in x-, y-, and g-directions, respec- 

tively, inches 

x, y, z rectangular coordinate axes 
a coefficient of thermal expansion, inches per inch 

per degree Fahrenheit 

;S angular rotation used in correction cycle, radians 

<y shear strain, radians 

S, A increment 
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e . normal strain, Inches per inch 

6 angular rotation about s-axis, radians 

X rotation of special set of axes, degrees or radians 

p normal distance to panel on i grid line, positive in 

positive ^-direction, inches 

er normal stress, positive for tensile stress, psi 

t shear stress, positive in direction of associated co- 

ordinate axis when tensile stress on cross section 
is in positive rr-direetion, psi 

4> angle between normal line r and 2 -axis, degrees or 

radians 

angle between normal line p and y- axis, degrees. or 
radians 

Subscripts : 
i, j, k grid system 

x. y, z . coordinate axes 
v, w, 6 bay displacements 

0 initial value 

n cycle of iteration 

A prime refers to the principal shear axes and two primes 
refer to the principal inertia axes. A bar over a symbol 
indicates an average value at the center of a bay. 

NOTATION' 

The notation employed is illustrated in figure 6. The 
system adopted for designating parts of the s tincture is as 
follows:- 

Bulkheads divide the length of the structure into a number 
of bays. The subscript i is used to designate a given bulk- 
head or the bay between the i— 1 and ith. bulkheads. 

The stringers and panels in a given cross section form the 
basis of a grid work which can be used to designate these 
elements. These grid lines are not necessarily straight, 
parallel, or perpendicular but follow the panels. Those grid 
lines that are approximately parallel to the 2 -axis are desig- 
nated by the subscript j; those approximately parallel to the 
y-axis are designated by the subscript k. 

With this system, points and stringers can be uniquely 
located as follows: 

The point on the ith bulkhead at the intersection of the 



yth and kth grid lines is designated by the subscripts ij,k. 

The stringer in the ith bay at the intersection of the jth 
and kth grid lines is designated by the subscripts i,j,k. 

In order to locate a panel, the grid line on which it lies 
must be known. This notation consists of underl ining the 
appropriate subscript; for example: 

The panel in the ith bay on the jih grid line and between 
the k — 1 and kth grid lines is designated by the subscripts 
i,j,k. 

The panel in the ith bay on the kth grid line and between 
the j — 1 and jth grid lines is designated bv the subscripts 
i,j,k. 

The grid lines and bulkheads are numbered such that the 
numbers increase in the directions of the positive coordinate 
axes. 

STBESS-STBAtM RELATIONSHIPS 


The shear strain in a given panel is constant and is 
defined by the following relationships which depend upon 
the location of the panel: 

=2 ~b f t *-!,*) + 


Aiq , A w t . , A 6 t 

JT COS jr sm j. ^ * 


(Ala) 




2hj t U ‘‘' 1 ’ k ~ 1 

AV{ . , A w t . . . A0, 

-j- sm i — cos Pit 


(Alb) 


When the shear strains are being- computed, the normal 
distances r. and p must be given their proper signs. 

The constant- shear stress produces a linearly varying 
strain in the stringer and its average value at the center of 
the bay is 


ti.i.t 


Vt ,},* — / 7 

L t \AE 



(A 2) 


Note that the thermal expansion is included in the relation- 
ship between stringer stress and strain. 


EQUEUBHrUM OF INDIVIDUAL STRINGERS 

If a half-bay length of stringer on each side of point (i, j, k ) 
is isolated, the force system of figure 7 is obtained, and the 
following equilibrium equation can be written: 


+ (rtL \ j _frg^ + f rtL^ _ 


+ 


\ 2 A,/.* V 2 \ 2 A r /.jt V 2 A^ r jt-{-i 


+i 


(rtL\ 




(it) +Cr) -(w) it; 

\ * A-t- l„/ r * V & A+1./-J-I,fc \ ^ /*+!,£* \ * A+l,J,*+l 

— P t. j.k= 0 (A3) 
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Substituting equations (Al) and (A3) into equation (A3), yields the following equation for the stringer displacement 
of point (i,j, k) in terms of the displacements of the adjacent bays and stringers: 


where 
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Equation (A4) involves no assumptions regarding equality 
of structural dimensions, temperatures, or elastic properties 
about point (i, j, k ) . If any element is missing, the associated 
stiffness goes to zero and tlie general equation is still 
applicable. Since AE and Gt always appear as products, 
the variation of elastic properties with temperature is 
equivalent to changes in the stringer areas and sheet thick- 
nesses of_the effective structure. Furthermore, the thermal- 
expansion terms appear in the same manner as axial loads 
applied to the stringers. Thus, if desired, tho effects of a 
' nonuniform temperature distribution can be determined by 
applying a set of equivalent loads to a new effective structure. 

BAY SHEAR AND TORQUE EQUILIBRIUM 

The equations for the bay displacements (t>/itf,0) can he 
obtained from equilibrium of the shear forces on the. bay 
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cross section 


[(jib cos — (rth sin ^)(./.*]=0 (A5a) 

it 

[(rtl sin <j>) { .j,k+(Tth cos $£,/,d=0 (A5b) 
Kjibr)i'j't — ]=0 (Aoc) 

Ik 


Substitution of equations (Al) for the shear stresses in equa- 
tions (A5) results in 


(«7»«)ih0i+CJ r tw)tAw { — ( Je,)t^6t=(F t ) c + ( u t.i.t-r^i-i,}.i) 


| ~(6t cos <f>\ / 

fOt cos <f\ 

Lv 2 Xw ' 

X 2 Jtj.t 




(A6a) 


K Ot sin <ft \ / Gt sin <f> \ 

2 /ui+ i,t \ 2 Jt.i.t 

(^X,, + r(^XJ (A6b > 

- % 

— — (t7su)jAiP £ -f-(J« s) £ A0 { =(Af c )f — 


where 


(J ) — sin <j> cos ^ sin ^ cos ^ J 

( j.j,= [(^^0..« + ( 2 T 2 ^X i ,] 

(AA '?f[(T0u,+(Tl t .] 

Equations (A6) can be simplified by elimina ting the cou- 
pling terms if the axes used in the computations are the prin- 
cipal shear axes of the cross section. These axes are defined 
such that 

0 (A7) 


The relationship between the location and orientation of 
points and panels in two systems of coordinates, arbitrary 
axes (x,y,z) and the principal shear axes (x',y',z?), is shown 
in figure 8 and given by the following equations : 



y'—(z—m r ) sin X'-f-f y—V) cos 

CA8a) 

sf ={z—m') cos X'—(y—l') sin 

“ (A8b) 

4>' —4> — X' 

(A9a) 


(A9b) 

r'=r+l' sin 4>—m' cos 4> 

(AlOa) 

p'=p—l' cos ip—m' sin f 

(AlOb) 

Then the location of the principal shear axes is 


tan 2 \' — T 2J ™ 

Wff WttU 

(Alla) 

7 / 

^ J T T 2 

• (Allb) 

m J T T 2 

C/ g|v MW t/ B o 

(Allc) 

and, with respect to these axes, 


cos 2 X'-\-J WiI1 sin 2 X'+2 J ra sin X' 

cos X' 

(Al2a) 

Jua'=Juv sin 2 X'+J ua cos 2 X'—%J clB sin X' 

cos X' 

(Al2b) 

Jtt',— Jlt~\~l' Jim — m'Jt, 

(Al2c) 

F V '=F X sin X'-\-F v cos X' 

(Al3a) 
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F/=F, cos \'-F, sin V (Al3b) 

M x '=M x +m'F a -l'F t (Al3c) 


When referred to the principal shear axes, the equations 
for the bay displacements become 


2 

/ Gt cos <j>' \ / Gtsin 4>’\ [ / Gtsm4/' \ "11 

\ 2 \ 2 \ 2 


(A 14a) 




/ Gtsin<t>' \ . / Gt cos / Gt cos i^ \ 1) 

\ 2 /<„/,* \ 2 \ 2 /*,/.*_!) 

’(Al4b) 


(Al7b) 

(S’ +(£’/„), t ,)i — (£iV/ J )f+ 

(Al7c) 

where 

(A£) { = (!?/»)«= 

(EI„) t =Z^(AEy\ }it 

(EQ v ) i =T l ^ l (AEz) i . i , l {EI ax ) t =^^{AEyz) t , lrk 

These equations can be simplified by elimination of tho 
coupling terms if the computations are referred to tho equiva- 
lent pidncipal inertia axes of the cross section. These axes 
are referred to as equivalent because the variation of modulus 
of elasticity over the cross section is taken into consideration. 
These axes (x",y",z") are defined such that 


BAY THRUST AND MOMENT EQUILIBRIUM 

The equations obtained from equations (A4) and (A14) 
are sufficient in themselves to define completely the dis- 
placements of the structure. However, if the equations are 
solved by iteration, it is helpful to employ a periodic cor- 
rection cycle based on the gross equilibrium of axial loads 
in the cross section . .... . 


(^h-S2pCP)<.,.*=o‘ 

(Al5a) 

(rf,)«-SS(^)ueO 

t k 

(Al5b) 

( .5?*) H - S X (Ey ) <. i. * = o 

k 

(Al5c) 

It can be shown that equations (Al,5) are 

satisfied by the 


solution of equations (A4) and (A14); however, they are not 
likely to be satisfied by the displacement values obtained 
from any given cycle of iteration. In reference 2 it was 
demonstrated that convergence of the iterative process can 
be speeded if the displacement values are periodically cor- 
rected so that the stringer displacements satisfy equations 
(A15). 

The corrections applied to the stringer displacements are 
a planar distribution over the cross section and are deter- 
mined as follows: 

(«<. t. k)«+i = ( u ( , *)„ ■ + A« ti j, k (A1 6) 

where 

Alt < J t *=£W7 + /3 f ,< 3/j,i+£y,J Z},k 

Substituting equation (A16) into,(Al5) yields 
C AE) { du t -r{EQ ! ) i Q y )i (3 y , ( = ( LF X ) ( + 

(Al7a) 


E0/'=EQ v "=EI t /'= 0 

and then the following relationships are applicable: 
y"=(z— m") sin \" + (y—l") cos X" 
z"—(z—m") cos sin X" 


2 {EI„-AEl"m") 

0 EI n -AEm ni )-(EI t ,~AEF *) 

M EQi 

~AE 


(A 18) 

(A 19a) 
(A 19b) 

(A20a) 

(A20b) 


m 


ff 


E Qp 
AE 


(A20c) 


EI„" ={EI„- AEm"*) cos* + {EI .-AEl"*) sin* X"- 

2 (EI vt -AEl"m") sin X" cos X" (A21a) 

EItt'— (JEI n —AEm"*) sin 3 +(EI,-AEl" 2 ) cos* X"+ 
2{EI a -AEl"m ") sin X" cos X" (A21b) 
F X "=F X (A22a) 

M,"=M t sin cos X" (A22b) 

M,''=^M x cos sin X" (A22c) 


A further simplification of the correction cycle can be made 
by eliminating the load and temperature terms on the right- 
hand side of equations (Al7). This elimination can bo 
accomplished by iterating the difference between the exact 
solution and one which satisfies statics (equations (A15)) 
but not. necessarily continuity. The iteration of differences 
has an additional advantage in that smaller numbers, and 
consequently less work, are required to obtain a solution. 

An examination of equations (A17) indicates that they can 
be satisfied by a planar distribution of strain corresponding 
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to the elementary analysis of reference 1. Then the initial 
values of stringer displacements u can be defined as follows: 

(eo-Zi)^ j r t (A23) 

where 

(e aL)t, j,t— (Stt/Ood - (3t.i /r )ayrx r + (3s,i ,r )a z i* ' 

and, with respect to the equivalent principal inertia axes, 

(SWf/ ° 0= (^) 1 [ (F/,)f+ 9? (ABaT)i - '■*] (A24a) 

(A24b) 

(k/'V = (^J( s «' , )<+SS ( ^^' ,) uJ 

(A24c) 

The corresponding values of the bay displacements are 
obtained from equations (Al4) and (A23). 

Then the correctiomcycle equations applicable to the 
iterated differences are as follows: 

hi" = — ( (A u i.j.k — (A25a) 

&/= i. t— Atii-i, } ' t )(AEy")i' 1 * (A25b) 

Pi.t"~ — (A25c) 

Equations (A25) provide corrections to the stringer dis- 
placements u only. These corrections remove, any unbal- 
anced moment or thrust on the cross section but add 
unbalanced shear forces which are removed by correcting the 


bay displacements (v,w,6). The corrected bay displacements 
are obtained from the corrected stringer displacements by 
application of equation (A14). These two operations con- 
stitute the complete correction cycle that brings the stringer 
loads into equilibrium with the external loads without chang- 
ing the shear stress in any panel. 

REFERENCES 

1. HeldenMs, Richard R.r The Effect of Nonuniform Temperature 

Distributions on the Stresses and Distortions of Stiffened-Shell 
Structures. NACA TN 2240, 1950. 

2. Dufaerg, John E.: A Numerical Procedure for the Stress Analysis 

of Stiffened Shells. Jour. Aero. ScL, voL 16, no. 8, Aug. 1949, 
pp. 451-462: 

3. Southwell, R. V.: Relaxation Methods in Engineering Science. 

A Treatise on Approximate Computations. The Clarendon 
Press (Oxford), 1940. 


TABLE L— STRINGER PROPERTIES 


® 


® 

© 

® 

© 

© 

<S> 

® 

8 

f 

£ 

/ 

T 

E 

aT 

A 

ae' 

AE. 

L 

AEcrT 





* 

•* 


©X® 

® 

® 

©X® 


20 

1 

203 

9. 55640X10* 

2,728.5X10-1 

1.11 

10.60760X10* 

530. 3SO 

23, 943 


20 

2 

158 

9.89078 

2; 066.0 

L 15 

IL 37440 

563,720 

23,500 


20 

3 

156 

9.80937 

2^025.2 

.74 

7.33293 

368,647 

14,851 


120 

4 

162 

9.8655$ 

2,120.6 

L0I 

9.96424 

493,212 

21, 130 


f2G 

1 

210 

9.45386 

2,910.7 

0.98 

9. 26478 

463,239 

26,967 


20 

2 

173 

9.79376 

2,271.3 

l. 05 

10.28345 

514.173 

23,357 

2 

20 

3 

165 

9. 84635 

2,101.6 

.60 

6.49859 

324,930 

14,04t 


120 

4 

1ST 

9.S3339 

2,189.0 

.89 

8.75172 

437,586 

19,158 


f2G 

1 

231 

9.35168 

3,080.0 

0.59 

5.51926 

275,963 

16,999 


J20 

2 

1ST 

9.6970S 

2,46X3 

.62 

6.01219 

300.610 

14,816 


l 20 

3 

1/6 

9.7S0S1 

2,298.8 

.38 

3. 71652 

185,826 

8, 544 


IS 

4 

172 

9.80044 

2,257.6 

.54 

5.29224 

2S4.6I2 

11,943 


f20 

1 

244 

9.24233 

3,284.4 

0.46 

4.25147 

212,574 

13,878 


I20 

2 

202 

9.5S691 

2,672.7 

.52 

4.98520 

249,260 

13,324 

4 

120 

3 

185 

9.7X125 

i,m j 

.30 

2.91338 

145,669 

7,099 


Lao 

4 

177 

9.70674 

2JS28J3 

.42 

4.10203 

205,102 

9,513 


*.£-*(10.5 — G.QOI/iT'T— 0.0030 LSI T 1 ) XIO" 
•♦a r- ( L2JS2 r+0. 00152 T’) X I (M 


TABLE II.— PANEL PROPERTIES 
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TABLE III.— PRINCIPAL SHEAR AXES AND BAY DISPLACEMENT EQUATIONS 
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TABLE IV.— PRINCIPAL INERTIA AXES, CORRECTION CYCLE, AND INITIAL STRINGER DISPLACEMENTS 
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TABLE V.— STRINGER-DISPLACEMENT EQUATIONS 

(a) Load problem * 
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TABLE V— STRINGER-DISPLACEMENT EQUATIONS— Concluded 

(b) Temperature problem * 
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844763 

0 

0.464713 


0 



0.285389 

0.041300 

0.083653 

a! 041306 



1 2 

1M 

.068705 

.360812 


.102630 



.223512 

. 034S26 

.001676 

.OOOSGS 

.008304 


1 3 

727,614 


.247336 

.099742 

.148722 



.150584 

.049576 

.000252 

.000345 

^^■TnrTTTrVl 


l 4 

810,844 


.450163 

0 

.133991 



.281854 

0 

—.080504 

-.044487 

.008892 j 


1 

557,833 

0 

0.432153 


0 

0.124224 

0 ‘ 

0.319400 

a 061671 

0.062553 

0. 001671 



2 

690,387 

.060543 

.332994 , 


^^■B7TTrr^H 

.103161 

.049830 

.239938 

.051270 

.001338 

.001446 



3 

474,609 

.073475 

.240067 


.150062 

.151478 

.074588 

.160802 

.075475 

.000529 

.000887 

.003045 


4 

541,607 

.066602 

.421966 


.132740 

0 

.066138 

.312553 

0 

-.006602 

-.066138 



1 

246,976 


0.721414 

0.139293 

0 

0.139293 


0 

0 

a 139293 




2 

319.002 


.562455 

.1X0950 

.107822 

.110950 


0 

0 

■00312S 

0 

JMfltiTsi n?« i 


3 


.163217 

.343252 

.165158 

.163216 

.165157 


0 

0 

.001941 





240,923 

.148082 

.702635 

0 

.148682 

0 

0 

0 

0 

-.148682 

0 

HHiil 




j* 


JttH.fK'll • 


AkA 


4- (AEaT)t tf . 1 — G4_EotT)<« rfj f 


\ 44 / {HJH4 

Obedc: @+@+@+@+®+@+®+@=l 


foti 

+ (' 


\ (go: 

A«J., 4 V 46 . 
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TABLE VI.— MATRIX OF COEFFICIENTS 

(a) Load problem 



Mm 

Min 

Mm 

11141 

Uni 

u«l 


0 

0.081708 

0 

0 

A 146946 

0.034293 


.089978 

0 

.148072 

0 

.060142 

.143783 


0 

.083376 

0 

.100257 

0 

.033617 

H 

0 

0 

.148138 

0 

0 

0 

u*u 


0.040592 

0 

0 

0 

0.081101 

*ai 


.107500 

.073689 

0 

.089839 

0 

«ai 

HM 1 

.04158! 

.183143 

.049984 

0 

.083138 



0 

.073856 

.124043 

0 

0 




um 


0 

.099491 
. 247536 
.099742 


0 

.148722 

0 

. 146tfl8 


0 

. 049231 
.166584 
.049576 


0 

0 

0 

0 


0 

0 

0 

0 


AOS 

w 

A6 t ' 

Adi' 


0. 657179 

0. 127908 

0.000600 

-0.869324 

-652002 

.128254 

—.000931 

-.863065 

0 

o 

o 

o 

0 

0 

0 

0 


0 



tin i 


0 

0 

.006826 

.165840 


0 

0 

.100041 

0 


0 

0 

.033215 

.127330 


0 

0 

0 

0 


0 

065132 

-.041655 

0 


0 

— L 153883 
-.662239 
0 


0 



0. 170292 
.045047 
0 
0 

a 010233 
.162371 
.041302 
0 

0 

.07347! 

.240057 

.073004 

0 

0 

.049980 

.191307 

0 

a 079899 

0 

0 

.083460 

0 

.150062 

0 

0 

.062119 

0 

.099012 

0 

0 

.151478 

0 


0.091092 

.044412 

0 

0 

0.039660 

.103648 

.040817 

0 

0 

.074588 

.156802 

.075475 

0 

0 

.049032 

.110031 

0 4 

0 

0 

0 

0 

0 

0 

0 

-.070462 

-.051331 

0 

-.002529 

-.069463 

-.050636 

0 

-.002258 


0 

0 

0 

0 

0 

0 

0 

0 

.053002 

.121609 

-.002644 

-.840016 

.048107 

.122248 

-.004436 

-.839620 





Mill 

0 

0 

0 

0 

0 

-.626850 

-.625728 

0 

0 

-.044014 

| -.043660 


UllS 

a 198993 

a 084548 

0 

0 

0 

0 

— 0.750708 

-A 754346 

0 

0 

j 0.036210 

' 0.036220 

tlHl 

.097020 

.225184 


0 

0 

0 

-.017604 

-.616690 

0 

0 

.007546 

.007648 

tem 

0 

.087001 

.343232 

.108289 

0 

0 

0 

0 

0 

0 

-.000090 

-.000142 

Mjfi 

0 

0 

.105187 

.240070 

0 

0 

—.825728 

-.829902 

0 

0 

-.043660 

-.043716 


«lll 

0 

0.084348 

0 

0 

0 

Mu* 

.097020 

0 

.163216 

0 

0 


0 

.087001 

0 

.108289 

0 


0 

0 

.165158 

0 

0 

AWi' 

0 

0 

mmm 

0 

0 

Atft' 

0 

0 


0 

0 

Attn' 

0 

0 


0 

0 

AWi' 

151742 

*>-. 107930 


— . 135638 

0 

iW 

0 

0 

0 

0 

0 

AW 

0 

0 

0 

i 0 

0 

Adj' 

0 

0 

0 

i o 

0 

AW 

L 415914 

.260572 

, -.009705 

! -1.831930 

0 


0.026848 






































































































































































STRESS ANALYSIS OF STIFFENED-SHELL STRUCTURES 
TABLE VI.— MATRIX OF COEFFICIENTS—Conoluded 

(b) Temperature problem 


i 


Um 

Um 

If tu 

tfm 

tfm 

Uni 

tfm 

am 

ttax 

i 

i am 

0 

0.103730 

0 

0 

0.464713 

0.058705 

0 

0 

0 

0 

1 tfm 

.125231 

0 

.148072 

0 

..083663 

.350312 

.099491 

0 

0 

0. 

uni 

0 

.106157 

0 

. 134615 

0 

.070331 

.247836 

.089504 

0 

0 

UlU 

0 

0 

.148138 

0 

0 

0 

.099742 

.450163 

0 

0 

> un i 

0.345615 

0.051558 

0 

0 

0 

0.102630 

0 

0 

0.432153 

0. 050543 

| um 

.062214 

.270760 

.073669 

0 

.124959 

0 

.148722 

0 

.002553 

.332994 

vnt 

0 

.052815 

.183143 

.067113 

0 

.105208 

0 

.133991 

0 

.051885 

l U24t 

0 

0 

.073856 

.337547 

0 

0 

. 149318 

0 

0 

0 

1 

' oau 

0 

0 

0 

0 

0.285369 

0.033925 

0 

0 

0 

0.100373 

i Umt 

0 

0 

0 

0 

.041306 

.223512 

.049231 

0 

.124224 

0 

' uni 

0 

0 

0 

0 

0 

.034828 

.156584 

.044437 

0 

.103161 

j UM1 

0 

0 

0 

0 

0 

0 

.049576 

.281854 

0 

0 


0.319400 

.001871 

0 

0 


0.049830 

.286033 

.081276 

0 


0.083015 

.062214 

0 

0 

0.001119 

.001258 

0 

0 

-0.000120 

.000186 

0 

0 

-0.067502 

-.067113 

0 

0 

0 

.083653 

.041306 

0 

0 

.001076 

.000898 

0 

0 

.000252 

.000345 

0 

0 

-.089504 

—.044487 

0 

0 

0 

.062553 

.061671 

0 

0 

.001338 

.001446 

0.001740 

0.000104 

0.000818 

0.001824 

0.011800 

0.008304 

•0.007363 

0.006892 

0.005595 

a 002161 

am 

tfm 

tUtr 

nut 

tfm 

ttui 


C*i'_ 

A*s' 

Ar*' 


uui 
Uffi 
| Uoi 

' «m 


0.328816 

.007047 

-.000542 

-.338321 


.068602 0 
.421966 0 


0.721414 

.139293 


.066138 0 
.312563 0 


0.107822 

.562455 

.110950 

0 

0 

.163217 

.343252 

.165158 

0 

0 

.148682 

.702635 

0.107822 

a 

0 

a 

.163216 

0 

.110950 

0 

.148682 

0 

.165157 

0 


0.327720 

.007995 

.000849 

—.£36564 

0.327720 

.007995 

.000849 

-.336564 


0.326763 

.008681 


0.326763 

.008681 

.002350 

-.337794 


0.325706 

.009449 

.003966 

-.339140 

0.325706 

.009449 

.003986 

-.339140 
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REPOET 1043 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TABLE VII.— CORRECTION CYCLE 


(a) Load problem 



df.i" 

/W' 

't 

P,a" 


-0.054016 

0. 053248 




— 057021 

.069103. 



Win 

037341 

.037350 



KMl 

-.050740 

.060299 




ttju 


-0.053248 

0. 053741 


urn 


-.059103 

.058541 


un i 


037350 

. .036188 


Ut 41 


-.050299 

.... 051505 



Mill 



—0. 053741 

O 0523211 




—.058541 

. 061348 




—.036188 

. ^35S62 

UUl 



— . 051505 






U4U 




—0, 0523910 





- — . 061348 

1 




—.035852 

t£4tl 





.-.050480 

«<ri" 

5 

5 

■Is 

5 

6th cycle 

-85.25 

—62. 43 

-.85; li 

-97.77 

10 th cycle 

-.7000 

— L 198 

L 799 

-a. 156 

14th cycle 

.09770 

.1949 

.1069 

.2673 


(b) Temperature problem 



ta," 



a., >" 

. bU%" 

(3. l" 


a «. 

ttm 

ttm 

ttm 

1*141 

-0.270067 

-.289578 

180687 

-.253677 

0.029722 

.009506 

-.008325 

—.030952 

0. 286240 
.295514 
.186749 
.251497 

-0.029018 

-.009780 

.008413 

.030984 

• 5i r 














— 

— _ 



III! 



-a 266240 
-.295514. 
-.186749 
—*251497 

0.029618 
.009780 
-. 008413 
-.030984 . 

a 268705 
' .292703 
.180939 
.257653 

-0.029564 

-.009645 



— 



_ 

. 

- — 



irausl 



— 

win 

ttui 

ttm 

ttwi 





*07208706 

•-.-293703 

t-.180939 

•*257053 

a 029564 
.009645 
-.007983 
-.031226 

0.201603 

.306742 

.179202 

.262100 

-0.020354 

-.010188 

.008130 

.031303 











’■ • ■ 






ttm 

ttAll 

ttill 

ttm 


y- 




_ 

-a 261*95 
—.306742 
-.179262 
-.252400 

a 029254 
.010188 
-.008139 
-.031303 



' 

. 

- 

■ 

— 

— 




— — 

— 

*.»" 

w,. 

«,<" 


—14.239868 
-4.239858 
6. 700142 
16.700142 


iill 

7 1 w 


’ -14. 275390 
-4. 275390 
5.724610 
15. 724610 


-14. 224G72 
-4.224672 
5.775328 
16.775328 








" 


— 

— 

— 

6th cycle 
9th oyde 

-83 

0 

JL 951 
-.0199 

1 

—64 

0 

3.713 

.05887 

-94 
~T „ 

.07708 

-98 

-2 

6.004 

.00613 





TABTJO VITI.— SUOCJSHSIVE VALUES Oil’ T)IHPTiAf! KMTSNT 

(a) Load imiblom 



Initial 

values 

1st oyolo 

DllToronoo 

2d oyolo 

3(1 oyole 

4th cycle 

6th oyole 
(oorrootlon) 

0 th oyole 

7th oyolo 

8th oyolo 

0th oyole 

10th oyolo 
(oorrootlon) 

nth 

oyolo 

12th 

oyolo 

13th 

oyolo 

14th oyolo 
(o irrcotlon) 

Totsl value 

Olieok oyolo 


-0.010205 

-0.017684 

-1280X1O" 4 

-1471 

-1346 

-2034 

-2210 

-2283 

-2318 

-mj 

-2341 

-2345 

—2348 

-2360 

-2361 

-2361 

—0.018640 

-0,018046 1 


-.010206 

-.010010 

-321 

-183 

-300 

-376 

-881 

-670 

-679 

—688 

-m 

-689 ' 

-690 

-500 

-690 


-.016886 

-.010886 


-.010208 

-.013036 

2069 

8002 

3407 

3610 

3334 

8808 

3428 

;m*u 

3447 

3443 

3447 

3449 

3449 


-. 612846 

-,01284ft 

«MI 

-.010206 

-.016304 

-0 

410 

021 

708 

629 

688 

618 

081 

637 

033 

037 

039 

<H0 

1140 

-.016066 

-.016084 


™Q. 020040 

-0,027070 

-1339 

-1710 

-2280 

-2643 

-2888 

-2060 

-8012 

-3030 

-3048 

—3054 

-3000 

-3062 

-3064 

-3068 

-a 029703 

-0.020703 


—.020040 

— , 026813 

-173, 

-10 

-183 

-267 

-600 

-609 

-012 

-018 

—621 

-027 

—629 

-629 

—636 

-920 

-.027260 

-.027209 


-.026640 

-.024000 

2031 

3760 

41fi2 

4812 

4000 

*100*1 

4136 

4184 

4108 

4167 

4101 

4163 

4L04 

4166 

-. 022476 

-.022475 

uu 1 

— . 026640 

-,020748 

-10H 

020 

007 


730 

831 

876 

806 

004 

m 

004 


007 


-.028732 

-.028732 


—0, 034430 

-0.036007 

-1177 

-1874 

-2101 

-23*16 

-2771 


-2046 

-2971 

-2083 

-2092 

-2999 


-8003 


-0,037432 

-0.037488 


-.034430 

-.034687 

-207 

-36 

-26 

-107 

-623 


-004 

-670 

-674 

—683 

-686 

-686 

—080 


-. 035116 

-.03611ft 


— 034480 

-. 032134 

2290 

3*180 

3908 

4100 

3734 

3837 

8884 

3908 

3914 

8906 

3910 

8912 

3013 

IKTiTfH 

-.030616 

-..03061ft 

MUI 

-.034480 

-.034423 

2 

738 

" 1144 

1324 

808 

090 

1061 

1074 

1083 

1074 

1080 


108*1 

1088 

-.038346 

-. 033344 


-0.030801 



-1232 

-1718 

-1940 

-2420 

-2808 

-2800 

-2684 

-2804 

-2606 

-2611 

-2013 

-2016 

-2014 

-0.030608 

-0.030606 


-.030801 


-83 

140 

07 

30 

-480 

-608 

-610 

-628 

—629 

-M0 

-641 

-641 

-642 

-641 

037432 

-.037432 


-.036801 


1066 

2086 

8802 

3478 

2980 

8094 

3133 

8140 

3150 

3146 

8160 

3161 

3152 

3163 

-.033738 

-.033738 

util 

-.030801 

-.030030 

-46 

832 

1104 

1383 

804 

077 

1021 

1041 

1049 

1038 

1042 


1048 

10*17 

-,036844 

-.036848 






HS1 


1676 


1710 

1720 


1736 

1737 

1737 



0.001176 

0, 00U74 








3700 


3802 

8817 


3846 

8847 

8849 



,110083 

.110060 






Ktrl 


3000 



4100 


4100 

4166 

4167 



.164026 

, 1.64026 

AK>(' 





1382 


3473 


8670 

8604 


8645 

3060 




.180310 

.100310 


ft ftOTCMfi 



-73 

-98 


-107 

-113 

-116 

-117 


-117 

-118 

-118 



-0. 000403 

-0. 000403 


— ! 000266 



-163 

-221 


—248 

-202 

-26ft 

-271 



-273 

-274 



-.000640 

—.000640 


—,000364 



-170 

-260 


-286 

-302 

-aio 

-314 


Mlti 

-317 

-317 


. 

-.000071 

-.000071 

m' 

-.000182 



-170 

-241 


-278 

-282 

-300 

-80S 


mi 

-800 

-307 



-.000400 

-.000400 


(b) Temperature problem 



Initial 

valuos 

1st oyole 

DUIorenoo 

2(1 oyole 

Sd oyole 

4th oyole 

8th oyolo 
(oorroo- 
tlon) 

6th oyolo 

7th oyolo 

8th oyolo 

oth oyole 
(oorroo- 
tlon) 

Total 

value 

Ohook oyole 

9m 

0,060286 

0.060710 

’ 464X10“* 

080 

734 

80S 

742 

780 

757 

7flfl 

-017 

786 

0.061012 

0.061011 

Will 

.046882 

.046079 

-603 

-723 

-810 

-856 

-807 

-014, 

-017 

-017 

.045666 

,046604 

HlJl 

,042907 

,042785 

-172 

—346 

-884 

-412 

-43-1 

-436 

-485 

-435 

—486 

.042472 

.042472 

«141 

.039232 

,080648 

317 

402 

522 

660 

664 

604 

602 

662 

502 

.030704 

.030706 

ttau 

0. 10*1043 

0.106679 

1030 

' 1288 

1644 

169*1 

1477 

1487 

1487 

1487 

1480 

0; 106120 

0.100120 

t*Mi 

,090486 

.096376 

-1110 

-1386 

-1694 

-1604 

-1714 

-1757 

-1750 

-1750 s 

-1750 

,004727 

.094727 

uni 

,088330 

.087856 

-474 

-723 

-H40 

-862 

-806 

-896 

-806 

-805 

ure 

-896 

.087436 

. 087436 

1*841 

,080174 

,080928 

761 

043 

1080 

1120 

1115 

1118 

1118 

1110 

.081203 

.081202 

1*311 

0.102000 

0.103640 

560 

1664 

1740 

1761 

1682 

1588 

1680 

1689 

1680 

0. 104570 

0. 10*1577 


.149633 

. 140052 

-681 

-1411 

-1682 

-1761 

-1868 

-1872 

-1872 

-1872 

-1874 

.147750 

.147760 

1*U1 

.130280 

,130048 

-238 

-831 

-885 

-800 

-003 

-062 

-061 

-001 

—068 

1106 

, 136328 

. 135323 

i*an 

.122034 

,123322 

388 

1008 

1107 

1204 

1103 

1195 

1100 

1100 

.124129 

. 124130 

1*411 

0,226499 

. 0.227888 

2389 

2086 

3054 

3080 

•' .'2882 

2866 

2867 

2857 

2864 

0.228383 

0.228868 


.206120 

.203755 

-2874 

-20J6 

-3088 

-3131 

-3253 

-3252 

-3262 

-3262 

-32M 

.202876 

.202875 

1*431 

■ ,186708 

.186693 

-1176 

-1701 

^1704 

-1724 

-1700 

-1780 

-1780 

-1780 

-1701 

.184077 

.184077 

1*W 

,107402 

’ .169166 

1763 

2148 

2160 

2164 

2144 

2146 

2147. 

2147 

2140 

.100648 

.100548 

Adi' 

Adi' 

Ana' 

Am' 

0.003676 

.011832 

.02150H 

.032718 



64 

144 

214 

330 

81 

184 

270 

870 


52 

143 

184 

288 

64 

147 

188 

266 

64 

147 

188 

266 




0.003726 

.011073 

.021696 

,032972 

0.003728 

.011077 

.021005 

.082072 


II' 

1 1 


I -i: 


■' *■! 
. ' j.d 


■!' 't 

''■] i*-‘ 



’ "Jil , 

!Jtr“ - i." 


r 'it 





% 


HI < ; I 


MU 
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TABLE IX.— STRESSES 


Load problem* 


Temperature problem** 


@ 

© 

© 

© 

© 

© 

@ 

© 

© 
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